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Abstract 

It is known that Heston's stochastic volatility model exhibits moment 
explosion, and that the critical moment s+ can be obtained by solving 
(numerically) a simple equation. This yields a leading order expansion 
for the implied volatility at large strikes: aBs{k,T)^T ~ — 1) x fc 

(Roger Lee's moment formula). Motivated by recent "tail- wing" refine- 
ments of this moment formula, we first derive a novel tail expansion for the 
Heston density, sharpening previous work of Dragulescu and Yakovenko 
[Quant. Finance 2, 6 (2002), 443-453], and then show the validity of a re- 
fined expansion of the type UBs{k, T)^T = {Pik^^'^ + /32 + ■ ■ ■)^, where all 
constants are explicitly known as functions of s+, the Heston model pa- 
rameters, spot vol and maturity T. In the case of the "zero-correlation" 
Heston model such an expansion was derived by Gulisashvili and Stein 
[Appl. Math. Optim. 61, 3 (2010), 287-315]. Our methods and resuhs 
may prove useful beyond the Heston model: the entire quantitative anal- 
ysis is based on affine principles: at no point do we need knowledge of the 
(explicit, but cumbersome) closed form expression of the Fourier trans- 
form of log St ( equivalent ly: Mellin transform of St)', what matters is 
that these transforms satisfy ordinary differential equations of Riccati 
type. Secondly, our analysis reveals a new parameter {"critical slope"), 
defined in a model free manner, which drives the second and higher order 
terms in tail- and implied volatility expansions. 

1 Introduction 



The Heston model [H] is one of the most popular stochastic volatility models 
used in the financial industry. Furthering its understanding, and in particular 
the understanding of its implied volatility surface, is of particular interest in the 
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light of the recent financial crisis: the volatility smile (underlying: SPX) did 
steepen after September 2008, then flattened again; it also steepened substan- 
tially after the flash crash in April of 2010 and has since flattened agair0. It is 
also worth recalling that the very existence of the volatility smile as we know it 
was triggered by the events of 1987. 

This general motivation is complemented by an everyday question in the 
financial industry: how to (smoothly) extrapolate the smile seen in the market 
(typically a stepping stone towards the robust construction of a local volatility 
surface). Theorem 3 below contributes precisely in this direction and we derive 
new expansions for the implied volatility in the Heston model. Recall that its 
dynamics under the forward measure are given by 

dSt = StVVtdWu 5o = l, 

dVt ^ {a + bVt)dt + cy/VtdZt, Vo = vq > 0, (1.1) 

where a > 0, 6 < 0, c > 0, and d{W,Z)t = pdt with p e [-1,1]. Observe 
that our choice 5*0 = 1, as well as zero drift, entails no loss of generality. As is 
well-known (cf. [U EJ [151 [22l [25] ) , the Heston model, as many other stochastic 
volatility models, exhibits moment explosion in the sense that 

T*(s) = sup{< > : E[Sl] < oo} 

is finite for s large enough. (Here and throughout the paper, E[-] denotes the 
risk-neutral expectation.) Differently put, for fixed maturity T there will be a 
(finite) critical moment 

s+ := sup {s > 1 : E[S^] < oo} . 

(In the Heston model, and many other affine stochastic volatility models, T* is 
explicitly known. The critical moment, for fixed T, is then found numerically 
from T*(s+) = T.) A model free result due to R. Lee, known as moment formula 
(cf. [1 [13]; see also [1 [H [H]), then yields 

limsupaBs{k,T)^T ^ ^{s+ - 1) x fc, (1.2) 

A;— f oo 

where k = log{K/So) denotes the log-strike, ass the Black-Scholes implied 
volatility, and 

•9{x) = 2-4:{\/x^+x-x) e [0,2]. 

We remark that, subject to some "regularity" of the moment blowup (fulfilled 
in all practical cases; cf. [2]), the limsup can be replaced by a genuine limit. 
Thus, the total implied variance aBs{k,T)'^T is asymptotically linear in k with 
slope ^'(s+). (Similar results apply in the small strike limit k — > — oo, but the 
focus of this paper is on fc — )■ oo.) 

Parametric forms of the implied volatility smile used in the industry respect 
this behavior; a widely used parametrization is the following. 

^From a private communication with a derivative trader at a major investment bank. 
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Example 1 (Gatheral's SVI parametrization [17]). For fixed T, a parametric 
form of aBs{k,T)^T is given by 



fc a + b 



(-m + fc) r + J (-m + k) + s 



= SVI(/c; a, b, r, m,s) 



An expansion for fc — > oo yields 

SVI(fc) = fc b (1 + r) + (a - bm (1 + r)) + 0(fc"^), 



+ 0{k-i), (1.3) 



and we see that SVI(A:) is asymptotically linear. Remark that this parametriza- 
tion is not ad-hoc but has been obtained fey a T — > oo analysis of the Heston 
smile; cf. fJW and fT?J. 

Our main results are the following two theorems. Remark 1151 in Section [3.31 
and formula (|4.1ip in Section 2] complement them by left-tail asymptotics. 

Theorem 2. For every fixed T > 0, the distribution density Dt of the stock 
price St in a correlated Heston model with p < satisfies the following asymp- 
totic formula: 

Dt{x) = Aix-^'e"^^^^ {log x)-^'^^'''''^ (l + 0((log.T)-i/2)) (14) 

as 2: —> 00 . The constants A3 and A2 are expressed explicitly in terms of critical 
moment s+ and critical slope 



dT*{s) 



as 



As = 



ds 



and 



s— s+ 



A2 = 2- 



(1.5) 
(1.6) 



An expression for Ai is presented in Remark \lS\ below. 



Theorem 3. Under the assumptions of Theorem\^ the Black-Scholes implied 
volatility admits the expansion 



as fc —>■ 00, where 
Pi 

/33 



TfT=[p,k^'^ + P2+P,'^ + 0{±,))' (1.7) 



M f_i 1 

v/2 

1 /I a 



V2 V4 



1 

VA^ 
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Remark 4. The restriction to p < is (mathematically) not essential, but 
allows to streamline the presentation. As is commonly noticed, this covers 
essentially all practical applications of the Heston model. We also note that, 
since {a + bVt) = —b{a/{—b) — Vt), it can be helpful to think of —b (resp. 
V = a/ (— &)) as the speed of mean-reversion (resp. mean-reversion level) of the 
Heston variance process. 

Let us draw attention to the main predecessors of this paper: Dragulescu- 
Yakovenko [5| apply a saddle point argument to deduce the leading order behav- 
ior of the density in the stationary variance regime; essentially Dt{x) ~ x~^^. 
Gulisashvili-Stein [20] study the "uncorrelated" Heston model {p — 0) and find 
the same functional form as in (|1.4|) and (|1.7|) . with (more involved) explicit 
expressions for Ai,f3i. (Their method relies on representing call prices as aver- 
age of Black-Scholes prices and does not apply when p ^ 0.) While it is easy 
to see that, in the case p — 0, our expressions for A3 agree, it is checked in 
Appendix H (for the reader's peace of mind) that our A2 = 2-^^-|p=o coincides 
with their expression for ^2- In Appendix HI we present a numerical example 
that shows the accuracy of our asymptotic formula for the density, and of the 
resulting implied volatility expansion. 

An interesting feature of our approach, somewhat in contrast to most ana- 
lytic treatments of the Heston model0 is that our entire quantitative analysis 
is based on afhnc principles; at no point do we need knowledge of the (explicit, 
but cumbersome) closed form expression of the Fourier transform of log St or, 
equivalently, the Mellin transform of St- (With one inconsequential exception, 
namely a simplification of the formula for the constant factor Ai.) Instead, we 
are able to extract all the necessary information on the transform by analyzing 
the corresponding Riccati equations near criticality, using higher order Euler es- 
timateslfl In conjunction with a classical saddle point computation we then "im- 
plement" the Tauberian principle that the precise behavior of the transformed 
function near the singularity (the leading order of which is exactly described by 
the critical slope!) contains all the asymptotic information about the original 
function. At this heuristic level, we would expect that the critical slope a, as 
defined in (jl.Sp . is the key quantity that drives the second and higher order 
terms in tail- and implied volatility expansions of general stochastic volatility 
models (even in presence of jumps). Back to a rigorous level, it appears that 
the key ingredients of our analysis are applicable to general affine stochastic 
volatility models (cf. [22J), and we will take up on this in future work. 

The explicit constants A^, /3i for z = 1, 2, 3 in the above theorem are clearly 
tied to the Heston model itself. In fact, it is the explicit nature of how these 
constants depend on the Heston parameters (a, &, c, p), as well as spot vol vq 
and maturity T, that furthers our understanding. Let us be explicit. It follows 
from equation p.4|) below that s+ = s+{b, c, p,T) does not depend on a,vo 
(equivalently: does not depend on v, vq); furthermore s+(T) — ^ 54.(00) £ (1, 00) 
as T — 7- 00. Moreover, the critical slope is explicitly computable: a/T will 

2 Exceptions include [T0l[22l . 

■^See 1161 for more information on the power of Euler estimates. 
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be seen to be an explicit fraction involving only b, c, p and s+ but not a, vq 
(equivalently: v,vo)- We see furthermore that l/cr — (T/a)/T — 0{1/T) as 
T — >• oo. As a consequence of all this, we see that changes in spot vol are 
second order effects: j3i does not depend on ^/vq, whereas /32 depends linearly on 
it. Practically put, we see that increasing spot vol allows to up-shift the smile 
(intuitively obvious!) but does not affect its slopes at the extremes. We also 
note that changes in v are not seen until looking at (3^. No such information 
could be extracted from (|1.2I) and previous works. 

Another application concerns the design of parametrizations of the implied 
volatility: the SVI expansion (|1.3I) is not compatible with the correct expan- 
sion (|1.7p : the latter has a constant term, /32, which is not present in (|1.3p . (We 
are grateful to J. Gatheral for pointing this out to us.) The solution to this 
apparent contradiction (recall that SVI was obtained by a T — oo analysis of 
the Heston smile) is simply that h A2 = 0{a~'^/'^) = ©(T-^/^) 0. In fact, 
this suggests that SVI type parametrizations could well benefit from additional 
terms corresponding to such a /32-term; essentially accounting for the fact that 
T^oo. 

2 Moment explosion in the Heston model 

2.1 Heston model as an afRne model and moment explo- 
sion 

Consider the correlated Heston model given by , and set Xt — log St ■ From 
basic principles of affine diffusions (see, e.g., [25]) we know that 

logi?[e^^*] =0(s,t) + z;oV^(s,i), (2.1) 

where the functions 4> and "0 satisfy the following Riccati equations: 

= i^(s,7A), </)(0)=0, (2.2) 
V' = i?(s,0), 7A(0) =0, (2.3) 

with F{s, v) = av and i?(s, v) = i(s^ — s) + ^c^u^ + bv + spcv. In (|2.3p . and ijj 
are the partial derivatives with respect to t of the functions and ip, respectively. 
Our goal in Section [2] is to identify the smallest singularity, s = s+, of (|2.ip . and 
to analyze the asymptotic behavior of (|2.ip in its vicinity. The estimates found 
will be put to use in Section |3l where we perform the asymptotic inversion of 
the Mellin transform £'[e'"~^-''^*] of the Heston model. 

Remark 5. The symbol s denotes a real parameter. The Riccati ODEs in p.2p 
and (|2.3p are also valid when s is replaced by a complex parameter s + iy. 

Given s > 1, define the explosion time for the moment of order s by 
T*(s) ^ sup {t>0: Eie"^'] < 00} . 
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An elementary computation gives 

2c^ min R (s,r]) ^ - \{spc + bf - (s^ - s)] -A (s) . 

176(0,00] L J 

Let us also set x(s) = spc + b. A typical situation in applications (a correlation 
parameter satisfying p < 0, and a non-zero mean reversion & < 0) implies that x 
is negative for s > 0. We thus assume in the sequel that 

x{s) < for aU s > 0. 
This assumption allows to use the following formula from |22, Theorem 4.2]: 
T*(s)-[ ifA(s)>0 

Remark 6. The integral in p.4p can be represented as follows: For A(s) < 0, 
we have 

T*(s) ^ — I arctan f^}^'^^ + tt] . (2.5) 

The derivative ^ ^ ^ 

dsT* = — ^is,v)dv 

can be computed explicitly. Indeed, from (j2.5p we get 

2pc{spc + b) - c^{2s - 1) 



2A(s) 

[c2(2s - 1) - 2pc{spc + b)] {spc + b) + 2pcA{s) 



A(s) [(spc + 6)2- A(s)] 



(2.6) 



2.2 Moment explosion 



For t > 0, let s+{t) > 1 be the (generalized) inverse of the (decreasing) func- 
tion r*(-), that is 

s+(i) = sup {s > 1 : E[e'^'] < 00} . 

Definition 7. Given T > 0, we call 

s+ := s+{T) = sup{s > 1 : E[S^] < 00} 

the "critical moment". The quantities 

a -.^ -dsT*\s^ >0 and K:^dlT*\,^ 

are called the "critical slope" and the "critical curvature" , respectively. Note 
that s+, a, and k depend on T. 
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Since T*(s+) T, formula implies that 

dT* ^ , i?i 



(2.7) 



where 



Ri = Tc2s+ (s+ - 1) [c^ (2s+ - 1) - 2pc (s+pc + b)] 
- 2 {s+pc + b) P (2s+ - 1) - 2pc {s+pc + b)] 



Ape 



c^s+ {s+ - 1) - {s+pc + by 



and 



i?2 



2c2s+ (s+ - 1) c^s+ {s+ - 1) - (s+pc + by 



Remark 8. The critical moment s+ can (and in general: must) be obtained by 
a simple numerical root-finding procedure. 

Let s > 1. We know that T*{s) is the explosion time of ■(/'• On the other 
hand, using the Riccati ODE for tj}, we see that 



± _ R{s,i^) 



Since R{s,u)/u'^ — > 0^/2 as w -> oo, we obtain 



i^{s,t) 



1 



i{T*{s)-t) 



eiSttT*{s), 



(2.8) 



uniformly on bounded subintervals of [1, oo). Next fix T > 0. Then we have T = 
T*{s+) with s+ — s+{T). Since the function T* is continuously differentiable 
(and even C^) in s, we have 



T* {s) - T = T* {s) - T* {s+) 

= {s+ - s) {a + 0{s+ - s)) 
^ a {s+ ~ s) as s t s+, 

where a — —dsT*\sj_ is the critical slope. Hence 



(2.9) 



xP{s,T) 



{s+ — s) c^cr 



s+{T). 



(2.10) 



It follows from (|2.8p and (|2.10p that (f>{s,t) — aip{s,'d)d'd has a logarithmic 
blowup: 

cj,{s,t) ^-^ log {T*{s)-t) asttr*(s); 



or 



(^,r) 



2a 



log((s* 



as s t s+ = s+iyr). 



The following lemma refines these asymptotic results. 
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Lemma 9. For every T > and for s 'f 5+ = s_|_(T), the following formulas 
hold: 

^{s, T) = - - ^ + 0(,s+ - (2.11) 

2a, 1 2a, T 

</>(s,T) = —log + — log- 



+ ay^ (^^(g+, z9) - „ ^) j + - (2-12) 

Proof. The idea is to use (second order) Euler estimates for the Riccati ODEs 
near criticality; this yields the limiting behavior of V'(s, t) and (/)(s, t) ast^ T*{s), 
and we complete the proof using (|2.9p . More precisely, let us introduce time- 
to-criticality r = T*{s) — t, and set V'(s,t) = 'ipis,T*{s) — r). Observe that 
l/tp{s,0) = and 

where W{s, u) = ^ + {b + spc)u + ^-^-^u^. A higher order Euler scheme for this 
ODE yields 

(1/Vi)(s, r) = (1/?/')(s, 0) + W(s, 0)r + W(s, 0)W(s, 0)t^/2 + o{t^) 

2 

as T — > and s stays in a bounded interval. Since W{s, 0) = ^ and W'{s, 0) = 
6 + spc, we obtain 

c2 / b + spc 2^^"^ 
= -r(^l-^^r + 0(r 

It follows that 

= ^-^ + 0(r) (2.13) 
as r = r*(s) - 1 1 0. Note that 
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Hence we obtain 

I f rr^ 2 b + s+pc K 

C fJ (s+ — s) C C O 

as s t s+ = s+(T). For the expansion of (/)(s, t) = a'0(s, we find 

/•* / 2 \ 2a /■* 1 



2a 1 2a 



+ aj^ + (2.14) 

To see the last equahty, note that the integrand of 

T'is) / 2 \ 

^(«' ^) - ^T^Twl = OiT* is) - t) 



c2(r*(s) -?9) 

has an expansion resulting from (j2.13l) , which may be integrated termwise [B' . It 
now suffices to use p.9p and (|2.14l) to see that, as s f s+ = s+(T), formula (|2.12l) 
holds. ■ 

Remark 10. It follows easily from the proof that Lemma |9] also holds as s 
tends to in the complex plane, provided that 5R(s) < s+. 

3 Mellin inversion via saddle point method 

Our proof of Theorem [5] proceeds by an asymptotic analysis of i?[e*^"^^-''^^], 
where u is complex. This is the Mellin transform of the density of St- As noted 
in Section 12.11 above, we can represent it in terms of the functions (j) ip 
appearing in the Riccati ODEs: 

log£;[e(""^)'^^] = 0(u- l,T) + voij{u-l,T). 

The density can be recovered using the Mellin inversion formula, that is 

-1 /"S+ioo 

where L = log a;, provided that s is in the fundamental strip, s e (s_(T), s+(T)). 



Remark 11. The integral in p.ip exists, since its integrand decays exponen- 
tially at ±ioo (see Lemma 1181 in Appendix I). Moreover, if u — 1 is imaginary, 
then the characteristic function of the random variable Xt ~ log (St) decays 
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exponentially. It follows that Xt (and therefore St) admits a smooth density. 
Since St is (a component) of a locally elliptic diffusion with smooth coefficients, 
this can also be seen employing classical stochastic or PDE methods (see [7] for 
some recent advances in this direction). 

We will deduce the asymptotics of p.ip by the saddle point (or steepest 
descent) method [HI [12]. The main idea is to deform the contour of integration 
into a path of steepest descent from a saddle point of the integrand. In cases 
where the method can be applied successfully, the saddle becomes steeper and 
more pronounced as the parameter (x in our case) increases. We then replace 
the integrand with a local expansion around the saddle point. The resulting 
integral, taken over a small part of the contour containing the saddle point, is 
easy to evaluate asymptotically. Finally, it suffices to show that the tails of 
the original integral are negligible, in order to establish the asymptotics of the 
original integral. Our treatment bears similarities to Taylor expansions studied 
by Wright and to the saddle point analysis of certain Lindelof integrals [TT] . 
The type of the pertinent singularity (exponential of a pole) is the same in all 
cases. 



3.1 Finding the saddle point 

A (real) saddle point of the integrand in formula p.ip can be found by equating 
its derivative to zero. Since it usually suffices to calculate an approximate saddle 
point, we note that Lemma[9]and Remark 1101 implv the following expansion, as 
u ^ u* -.^ s+ + 1 = A3 with 5J(u) < u*: 



where we put /S^ — 2vQ/c'a and 



/32 



— u 



2a, 1 



+ T + Oiu* -u), (3.2) 



r = -Vo 



b + s+pc 



2a, T 



^log^+ay^ (^^(.„.)--,^jd.. 

(3.3) 

Retaining only the dominant term of p.2|) , we get the approximate saddle point 
equation: 

/?2 ■ 



or equivalently. 



X exp 



-L 



0, 



/32 



0. 



[u* — uY 

The solution to the previous equation, 

u = u(x) :— u* — j3L-'/', 
is the approximate saddle point of the integrand. 
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3.2 Local expansion around the saddle point 

Our next goal is to expand the function 0(u — 1,T) + VQ'ip[u — \,T) at the 
point u — u. Put u = u + iy, and recall that we use the following notation: a = 
—dsT*\g^ and L = log a;. Since the (approximate) saddle point u approaches u* 
as i — ?> cxD , we may find the expansion of the integrand using p.2p . To make the 
expansion valid uniformly w.r.t. the new integration parameter y, we confine y 
to the following small interval: 



(3.4) 



The choice of the upper bound on a in p.4p will be clear from the tail estimates 
obtained in Appendix I. Since u* ~ u — jSL'^^^'^ — iy, we have 



1 



u 

= r^L^'^il + ir^L^/^y - r^Ly^ + 0(L3/2-3")) 



It follows that 



(3.5) 



log- 



log 



= -l0gL-l0g/3 + 0(il/2-a) 

Next, plugging the previous expansions, with u = ii + iy, into p.2p . we obtain 
the following asymptotic formula: 

(j){u — 1 + iy,T) + voip{u — I + it, T) 

= + ^ log L + ^Ly - ^^L^/^y^ _ ^ log /? + T + 0{L^-^^). (3.6) 



3.3 Saddle point approximation of the density 

For the sake of simplicity, we will first obtain formula (jl.4p with a weaker error 
estimate 0((log a;)"^/^"^^), where e > is arbitrary. Then it will be explained 
how to get the stronger estimate 0((log x)~^/^). 

We shift the contour in the Mellin inversion formula p.ip through the saddle 
point -u, so that 



Dt{x) 



1 



u~ioo 
oo 



27r 



^-uL+,p{u-l,T)+voi>(u-l,T)j^^ ^3j-j 
Q-iyL+4>{,u+iy-l,T)+voil>{u+iy-l,T) fj^y ^3 



The term 



11 



will yield the leading-order decay in (|1.4p : its exponent corresponds to the lo- 
cation of the dominating singularity of the Mellin transform. The lower order 
factors are dictated by the type of the singularity at u = u* , to be unveiled in 
what follows. 

The "tail" of the last integral in (13. 8p . corresponding to \y\ > L"", can be 
estimated using Lemma EOl fsee Appendix I). Therefore, 



£)y^2;)=x~"— 1^ ^-iyL+tj}(u+iy-l,T)+VQii(u+iy-l,T) fj^y 



2tt 



X exp 



- /3-ii3/2-2" + O(logL)) 



Next, using p.6p and the equality x "exp(/3i^/2) — x ^ exp(2/3L^/2), we 
obtain 

X (1 + 0(i2-3")) + ^'A, (^p^in _ ^-1^3/2-2a _^ Q^j^g j _ (g g) 

Evaluating the Gaussian integral, we get 

exp(-/3-iL3/2y2^)^^ ^ ^1/2^^3/4 / exp(-w2)du; 

-L-° J_^-l/2i3/4-a 
/•CO 

^1/2^-3/4 / exp(-u;^)d«; = V^/3i/2^-3/4^ (g^^Q) 



Here we use the fact that the tails of the Gaussian integral are exponentially 
small in L. Taking into account p.9p and p.lOp . we can compare the main part 
of the asymptotic expansion and the two error terms: 

const X x'^'^L"/"^'^/^ exp(2^Li/2) (main part) 

^-A3^a/c^-3/4 exp(2/3ii/2) Q(^2-3a) (^j.^.^^. f^.^^^ j^^^j expansion) 

a;"^^ exp(2^L^/2 - ^-1^^/2-2" ^_ Q^i^g (gj-j-^j. fj.^^^ ^^il estimate) 

Since 2— 3q; < 0, the expression on the second line is asymptotically smaller than 
the main part. In addition, since 3/2 — 2q; > 0, the quantity exp(— /3~^L'^/2~2") 
decays faster than any power of L. This shows that the expression on the 
third line is negligible in comparison with the error term in the local expansion. 
Hence, it suffices to keep only the error term resulting from the local expansion. 
As a result, the error term in the asymptotic formula for Dt is ©(L^"^") = 
0(i-i/4+e), (Take a close to |.) More precisely, using p.9p and p.lOp . we get 
the following formula: 



Dt[x) 



exp(r) r^„i/2^2a/c^ 



a;-(«+ + l)e2/3i''"2,-3/4+Q/c= 



27r 

X (l + 0(L-i/4+=)). (3.11) 
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It follows from p. lip that formula (|1.4p , with a weaker error estimate, holds 
for the correlated Heston model of our interest. 



Remark 12. The integral on the right-hand side of p.3p can be easily calcu- 
lated from the closed form expression [SI HI] of tp. By (|3.1ip , we thus obtain 
the explicit expression 

( fb + s+pc K \ aT 
X exp -vo 5 h -ir-^ 5-(& + cps+) 



X 



' 2v/b2 ^ 2&cps+ + c2s+(l - (1 - \ ° ' 

c2 s+ (s+ - 1) sinh \ 26cps+ c2s+(l - (1 - p2)s^) j 



for the constant factor in (11.41 



Our next goal is to show how to obtain the relative error 0((loga;) ^^2) jn 
formula ()1.4p . Taking two more terms in the expansion p.Sp of 1/(m* — u), we 
get 

= + ,^-2^y _ ^-3^3/2y2 _ j^-4^2^3 _^ ^-5^5/2^4 _^ ^(X^-S"). 

Expanding the logarithm, we obtain 

log - iog(r + ir^L^'^y - r^Ly^ + ©(i^/'-^"))) 



u — u 



= i logL - log/3 -f ^i^/'y - \r'Ly^ + 0{L^/^-'-). 



We insert these two expansions into p.2p to obtain a refined expansion of the 
integrand: 

x-^-'y exp (0(w - 1 + iy, T) + wo^/;(u - 1 + it, T)) 
- X--' exp (^2/3^1/2 ^ log L - r'L^^V " J log /? + 

(l + ciL^y^ + C2L'^V + cL^'^y + C4 V + ^5^-^^' + 0(i-^/^+^; 

(3.12) 

for some constants ci, . . . ,05. Note that the terms with ci and C2 come from 
{u* — u)^^ , those involving C3 and C4 from log(u* — u)~^ , and the one with C5 
from u* — u. (To be precise, we have used that the 0()-term in (|3.2p is of the 
form c{u* — u)-\-0{{u* — u)2), as is easily seen by a third order Taylor expansion 
along the lines of Section 
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We will next reason as in the proof of the weaker error estimate. The main 
term and the error term from the tail estimate remain the same. The error term 
from the local expansion can be obtained as follows: Integrate the functions on 
both sides of formula (|3.12l) and take into account that 

^ 2/3gxp(-/3-'i'/'y')dy = y^ yexp(-/3-ii3/2y2^dy = 0. 

The two integrals resulting from the and y^-terms in p.l2p are easily calcu- 
lated; they yield a relative contribution of i^^/^, which merges with the term 
CsL"^/^. Hence we see that the absolute error term from the local expansion is 

^_A3^a/c^-3/4g^p(2/3Ll/2) X 0(L-l/2)_ 

This completes the proof of Theorem 

Remark 13. Note that the preceding argument can be extended by taking more 
terms in the local expansion of the integrand. A full asymptotic expansion in 
descending powers of i = log x can thus be obtained, which replaces the error 
term (1 + 0{{logx)-^/'^)) in ([13) by 

1 + Ci(loga;)-i/2 _^ C'2(logx)-3/4 + . . . + 0((log x)-™/^) 

with some constants Ck and arbitrarily large m. This is a typical feature of the 
saddle point method (see [12|, Section VIII. 3). 

Remark 14. By a standard result on integrating functions of regular varia- 
tion [SI Proposition 1.5.10], formula (|1.4p yields the estimate 

P[St > x] = -^x'^'+^e^'^^ilogxy'^^^+^^^'h+Oiilogx)-^^^)), 

as X 00, for the tail of the distribution of St- Note that the main fac- 
tor 3;"'^^+"'^ has been obtained by Dragulescu and Yakovenko [HI Section 6]. 

Remark 15. We briefly discuss the behavior of the Heston density Dt{x) near 
zero. Define the lower critical moment by 

s_ := inf {s < : E[S'j] < 00}, 

and the corresponding slope and curvature by 

o-_ := dsT*\,_ > and k_ := d^T*\,_. 

As X 4, 0, the integrand in p.l[) has a saddle point that approaches the singularity 
s_ -I- 1 at a speed of (— \ogx)~^^^ . All steps of the subsequent analysis precisely 
parallel the case x —i' 00 treated above. The net result is 

Dt{x) = Bix^^'e^'^^^{~logx)''^^'~^^^{l + 0{{-logx)-^/^)) (3.13) 
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as X i 0, where 



S3--(.s- + l), ^2 = 2- 



(2«o) 



1/4-a/c^ 2a/c^-l/2 -a/c^-1/4 



X exp —Wo 



b + S-pc K 



aT 

2^ . 2" & + cps_ 



2^/62 + 26cps_ + c2s_(l - (1 - p2)s_) 



2a/c 



- (s_ - 1) sinh i y/b'^ + 2bcps^ + c'^s^{l - (1 - p2)s_) 




Figure 1: — log-Dy^(a;) with its asymptotic approxiraations, where Dl^^ is the 
density of log St- 

Remark 16. The density _D^^ of the log-spot price log 5*2- is given by 
Its asymptotics readily follow from (|1.4p and (13.13^ : 



log ;c 



a; — oo, 



and 



Figure [1] shows the numerical fit of these approximations, using a set 

a^vX, 6 = -A, c = 0.2928, uo = 0.0654, /? = -0.7571, 

V = 0.0707, A = 0.6067 (3.14) 



of typical market parameters 
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4 Call pricing functions and smile asymptotics 



Recall that our main result (Theorem [2]) is the following asymptotic formula for 
the stock price distribution density in a correlated Heston model with Sq = 1: 

Drix) = Aix-'^'e^''^^i\ogxy^+^ {1 + Oiilogx)-^)) (4.1) 

as a; — oo. In the present section we will characterize the asymptotic behavior 
of the call pricing function K H' C{K) in such a model, and then prove Theo- 
rem [3] The following formula is a generalization of a similar result obtained for 
uncorrelated Heston models in [T5] : 



^ ' (-A3 + l)(-A3 + 2) ^ ^ ' 

X (l + o((logif)-^)) (4.2) 

as if — oo. Formula (|4.2p follows from (|4.ip . Theorem 7.1 in 19 , and Remark 
6.1 in [19]. Note that A3 > 2. 

We will next use the tail-wing formulas obtained in [H] to study the asymp- 
totic behavior of the Black-Scholes implied volatility K 1— >■ aBs{K,T) in a cor- 
related Heston model in the case where the maturity T is fixed and the strike K 
approaches infinity or zero. The following statement was established in (19j . 
Section 7. Suppose that the stock price density Dt in a general stock price 
model satisfies the condition 

cix-^h{x) < Dt{x) < C2X~^h[x) (4.3) 

for all large x, where ^ > 2, /i is a slowly varying function, and ci and C2 
are positive constants. Then for every positive function Lp on (0, 00) with 
livux^oo fix) = 00, we have the following: 



v T / 111 
aBsiK,T)^ = . /log A" + log - - log log- 



V2 V ^ K^Dt{K) 2 " K^DriK) 



- f °s low) ' 2 low) + (^^^^^^"^ ^^""^ 

= Jlogif + log-— ^— — - iloglog— , ^ 



K-i+^h{K) 2 ^ K-i+^h{K) 



' log TT^rTTTTTlTT ~ ^ log log ■ 



K-i+^h{K) 2 ^ ^ K-^+-^h{K) 
0({\0gK))-"^^(K)) (4.4) 



as if — )■ 00. 
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A similar assertion holds for small values of the strike price (see [TO], Sec- 
tion 7). It can be formulated as follows: Suppose that the stock price density Dt 
is such that 



for all sufficiently small a: > 0, where 7 > — 1, /i is a slowly varying function, 
and Ci and C2 are positive constants. Let t be a positive function on (0, 00) 
with liniK^o t{K) = 00. Then 



as 7^ ^- 0. 

Remark 17. The asymptotic formulas in (|4.4I) and (|4.6p are equivalent to 
similar formulas with (p{K) — 1 and t{K) = 1, respectively. Indeed, if for some 
function / and all functions g, which tend to infinity, we have f{K) = 0{g{K)) 
as K 00, then f{K) = 0(1) as if — 00. This can be shown as follows. If the 
function / is not bounded near infinity, then there exists a sequence Kn t 00 
such that f{Kn) > 2" for all n > 1. Put g{Kn) = n, and define the function 
g by linear interpolation. Then g{K) — cx) as if — 00, but f{K) 7^ 0{g{K)) 
as K ^ 00. The proof for A' ^ is similar. The authors thank Roger Lee for 
bringing this simple fact to their attention. 

Now let us apply (|4.4p and (|4.6|) to the Heston model. It is easy to see 
from (14. ip that (|4.3p holds with ^ — ^3 and the slowly varying function 



cia;^/i(a;-i) < Dt{x) < C2X^h{x~^) 



(4.5) 




(4.6) 
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It follows from (14. 4p and Remark [T7] that 



aBs{K,T)^ 



3\. , „ 1. . 1 



'{A3^l)logK-A2^logK- loglogX--loglo: 



a 3\. . 1. . 1 



47 2 ^ K-^^+'^h{K) 

0((logif)-t) (4.7) 



as if — >■ oo. Next, using the mean value theorem, we see that it is possible to 
replace the term i log log j^^Ag+a^^j^^) under the square roots in formula (|4.7I) 
by the term ^ log log K. Therefore, 



'{A^-l)\ogK~A2^/^^- ( ^-ijloglogif 



'(A3-l)logi^-A2yiog¥- ( __-)loglogif 



a 1 

? ~ 4 

+ 0(^(log/s:)-5j (4.8) 

as iiT — >• oo. Since \/l — = 1 — i/i + O (/i^) as /i 0, formula (|4.8p implies 
that 

.Vt 



aBs{K,T)^ = i^^As - 1 - - 2j (4-9) 



^ ^ / 1 1 



2 V\/A3-2 VA3 - 1, 

as K 00. Next, using (|4.10l) . we obtain the expansion (|1.7p for the im- 
plied volatility k i~->- aBs{k,T), considered as a function of the forward- log- in- 
moneyness k — log K. Theorem |3] is thus proved. In the case where p — 0, 
formula (|1.7p was obtained in [2D] (see and [12] for more details) . Note that 
already the leading order term 



aBsik,T)VT^pik'^^, 



00, 



gives very good numerical approximation results. This term was obtained in [2]. 
As a "lim sup" -statement, based on Lee's moment formula, it appears already 
in p. 
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Let us denote by Wbs the Black-Scholcs implied total variance defined by 
WBs{k,T) ^ <JBs{k,T)^T. 
Then formula (11.71) implies the following expansion for Wbs- 

WBs{k,T) ^ (3lk + 2(3iP2k^/'^ + 2l3il33\ogk + 0{ip{k)) as fc ^ oo, 



where (32, Pa, and (f are the same as in (|1.7p . 

Similar reasoning can be used in the case where k — — oo. Put j — B3 and 

h{x) = e^^^^(loga;)'^/^'-3/4^ 

where B2 and B3 are defined in Remark[Tni In addition, fix a positive function ip 
on (0,00) with linix^ao ^(x) = 00. Then (j3.13p shows that all the conditions, 
under which formula (j4.6p holds, are satisfied. Next, using (j4.6p and simplifying, 
we obtain the following asymptotic formula for the implied volatility in the 
Heston model: 

aBsik, T)Vf = pr(-fc)^/^ + P. + + O (4.11) 

as fc — >■ — 00. The constants in ()4.1ip are given by 

Pi = 72(7^3 + 2- /b^TT 

B2 f 1 1 

P2 = 

P3 = 



V2 \Vb^ Vb^, 

1 /I a\ / 1 1 



^/2 V4 c^J\VWT2 VS3 + I, 
For the total implied variance, we have 

WBsik, T) = pli^k) + 2pip2(-fc)^/' + 2pip3 log(-fc) + 0{ip{^k)) 

as k ^ —00. 



Appendix I: Tail estimates 

It is known 0|2S] that all the singularities of the Mellin transform _B[e("~^^'^'] 
of the stock price density Dt in the Heston model are located on the real 
line. Therefore, the function u i-)- e'^("~^'^)+'''''^("~-'^^^) is analytic everywhere 
in the complex plane except the points of singularity on the real line. The next 
statement justifies the application of the Mellin inversion formula in p.8p . and 
will be useful in the tail estimate for the saddle point method. By symmetry, it 
clearly suffices to consider the upper tail > 0). 
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Lemma 18. Let T > and 1 < si < 5ft(s) < S2- Then the following estimate 
holds as 3(s) — > oo; 



^^{s,T)+voip{s,T) 



where the constant C > depends on T , si, S2, and vq. 

Proof. Let s — £^ + iy and suppose y > 0. We will first estimate the function ip. 
Recall that 

1 

ip = -{s'^ ~ s) + y + bip + sippc with ^(C, 0) = 0. 
Set V-" = / + *.9 ^nd 7 = — (fo + Cpc). Then 7 > 0, and we have 

f-l{e~y'~^) + "^{f- 9') - If, f{s, 0) = 0, 

Our goal is to show that there exists a positive continuously differentiable func- 
tion t ^ C{t) on [0,T] such that 

f{s,t)<~C{t)y, (4.12) 

where s — + iy, 1 < si < < S2, and y is large enough. We first observe 
that / satisfies the differential inequality 

i<\{e-y'-^) + '^f-if (4.13) 

<~\v^ + \f-lf (4.14) 
for y > yo, where t/o depends only on si and S2. Set 

1 

Then (|4.14p can be rewritten as follows: 

f{s,t)<V{y,f{s,t)) (4.15) 

where s — + iy. 

We will next find a function C(t), t G [0, T] with C(0) = 0, strictly positive 
for i > 0, and such that the function F{y,t) —C{t)y satisfies the differential 
inequality 

Viy,F)<F. (4.16) 

Let us first suppose that such a function C exists. Then it is clear that given 
s = £, + iy, the initial data F{y, 0) — f{s, 0) = match. Now we can use the 
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ODE comparison results and derive from (I4.15P and (|4.16p that (I4.12p holds, 
which implies the following estimate: 



- ^vof{s,T) ^ ~voC{T)'i{s) 



(4.17) 



for all s = £, + iy with y large enough and si < ^ < S2- 
We now look for the function C satisfying the equation 

c{t) = ~jc{t) + e, 

where 6* is a positive constant, and C(0) = 0. The solution of this equation is 
given by 

'6*7-1(1 - e-T*) if7>0, 
et if 7 = 0. 



It follows that for t € (0, T] 



< c(t) < re. 



Next, choosing 6* > for which — i + = —j, we obtain 

Viy,Fiy,t)) <~ly'+ '^T'e'y^+-/Cit)y 

= -^y' + [e - C{t)) y 

<-C{t)y^F{y,t). (4.18) 

In (I4.18p . y is large enough and depends only on 9, and hence on the model 
parameter c and on T. This shows that the function F satisfies the differential 
inequality in (|4.16l) . and it follows that estimates (|4.12p and (|4.17p hold. 
Finally, we note that 



3?(0(s,r)) -a^ f{s,t)<ay(^ C{t)dt^ ^ ~ayC{T) 

Therefore, for 5(s) large enough, 

^<^{s,T)+v,^(s,T) I < exp |- (aC{T) + voC{T)^ 5(s)} . 
The proof of Lemma [18] is thus completed. ■ 



Lemma 19. If B > is any constant, then the portion of the integral p.Tp 
where 5(u) > B is 0(a;-^3 exp(^Li/2)). (Recall that L = logx.j 

Proof. If _B > _B is a sufficiently large positive constant, then it easily follows 
from Lemma [TSl that 



/ _ g-«i+0(«-i)+»'oV(«-i)rfy < Cx-^^ exp(/3ii/2) / g-cy^^ 

Ju+iB Jb 

= O (x-"^--" exp(/3ii/2 
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(The integral is clearly 0(1).) Moreover, since the Mellin transform of Dt does 
not have singularities outside the real line (see [26]), we have 



U+iB 



(,-uL+<t,+voi>^^ 



U+iB 



0(e-"^) = O (x~^' exp(/3Li/2)) . 



This completes the proof of Lemma [191 ■ 

Lemma [T9l shows that the part of the tail integral where 5(m) > i? is asymp- 
totically much smaller than the central part. We will next estimate the whole 
tail integral. 

Lemma 20. The following estimate holds for the tail integral: 

pu+ioo 



-uL+tp+votp 



du 



u+iL~ 



X exp 



(2/3L1/2 _ l^-1^3/2-2a ^ 0{\0gL)) . 



Proof. We will prove that there exists a constant -B > such that the absolute 
value of the part of the tail integral where L~" < S(u) < B equals 



_ 4 

X exp 



(2/3X1/2 _ l^-1^3/2-2a ^ O(logi)) 



(4.19) 



It suffices to establish this statement, since Lemma [TOl shows that the absolute 
value of the integral from ii + iB to u + ioo is asymptotically smaller than 
the expression in ([4J9|) . (Indeed: Dividing ([4l9l) by x"^^ exp{l3L^/^) yields 
exp(^Li/2 + 0(i3/2-2")), ^^[^-^^ tends to infinity Note that 3/2 - 2a < 1/2 
by 

It follows from Lemma [Hi and Remark [TUl that for some constant 7 > 0, 



^<j,{u-l,T)+voijiu-l,T) ^ Q 



exp 



13' 



^.s - u 



7 log (A3 - u) 



as u tends to w* s_|- + 1 = ^3 inside the analyticity strip. More verbosely, 
there exists a constant C > such that for a sufficiently small number B > 
and for all u in the analyticity strip with |5(u)| < B and 3?(u) > u* — B, we 
have 

13' 



- U 



Hence 



U+iB 



u+iL 



^ — uL+(t>+vat(: 



du 



< Ci-* exp(pL^''') £ ^ 1^3 - (a + exp (sR (-j-_^-_^^ 

= Cx-^' exp (2^2.1/2 _ ^-i^3/2-2a _^ o(logL) 
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We have used that the factor {A^ — {ii + iy)\ grows only hke a power of L, 
since 

pL-^ = A3 - ti < 1^3 - (m + iy)\. 
Furthermore, the quantity 

\A3-{u + iy)J (A3-u)2 + y2 

decreases w.r.t. \y\. Therefore, the integral of (|4.20p can by estimated by 
the value of its integrand at L^" times the length of the integration path. The 
latter is absorbed into C, and the former is given by 

(This can also by obtained by plugging y — L^°' into the singular expansion p.Sp 
computed above.) Finally, we write the factor L^^"^ as exp(0(log L)). ■ 



Appendix II: Comparison of constants 

Since s+ is the order of the critical moment, it is not hard to see that if p — Q, 
then the constant A3 defined by A3 = s+ + 1 is the same as the constant A3 
in [20]- 

We will next show that for p ~ 0, the constant A2 defined in (|1.6p is the 
same as the corresponding constant in [50] . It follows from (|1.6p and from (|2.7p 
that the constant A2 used in the present paper for p = satisfies 

A^ = ^ (4.21) 

with 

_ {2s+ - 1) [Tc2,s+ (,s+ - 1) - 26] 



2s+ {s+ - 1) [c2,s+ (,s+ - 1) - 62] • 

We will next turn our attention to the constant A2 in |20) . Lemmas 6.6 and 
7.3 established in [5D] provide an explicit expression for this constant. First note 
that the quantity r ~ fLT\b\ ™- |20j and the quantity 5+ in the present paper 
are related by 

r = |[cV(s+-l)-62]^ (4.22) 

This follows from the formula for A3 in (jl.6l) and from Lemmas 6.6 and 7.3 
in [2ni- 

It was shown in [20] . Lemmas 6.5, 6.6, and 7.3 that the following formula 
holds: 

r3(8C + T)i 
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with 



B = 



'2T I Tvq sin r 



c y2c2|| 1(1 + iT|6|)cosr-rsin 



2^/2y/vo\/ r sin i 



(l + iT|6|) 



cosr — rsmr 




and 



2c2 V 



Hence, 

c^VTy'2s+ - 1 1(1 + ^T\b\) cosr - rsinrj^ V 

Here we use the formulas for A3 in (|1.6I) and in Lemma 7.3 in [20]. Since 

lul 

between the constant A2 in [5D] and s+: 



r cosr + iT|6| sinr — and formula (|4.22p holds, we get the following relation 



ry2i7^[ir|6| (l + iT|6|)+r2]' 



4yi}^^s+ (s+ - 1) [c2s+ (s+ - 1) - 62] 



c^'2s+ - 1 [Tc^s+ (s+ - 1) - 26] 5 

XlicrcforG 

,2 ^ 16^0^+ (g+ - 1) [c's+ (,s+ - 1) - 62] 
2 c2 (2s+ - 1) [Tc2s+ (s+ - 1) - 26] ■ ^ ■ ^ 

Next, comparing (j4.2ip and (|4.23p . we see that the constant A2 used in the 
present paper coincides with the corresponding constant in [50]. 



Appendix III: Numerical results 

To conclude we illustrate the accuracy of (|1.4p by a numerical example, and show 
plots of the corresponding smile approximations. We will use the parameter 
values p.l4|) . Note that (fOI implies that 

\ogDrix) ^^^^ 33 2124, (4.24) 



log a; 
\og{x^^ Dt{x)) 
A/Iogx 



A2 w 12.3533, (4.25) 
Ai w 2311.69, (4.26) 
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as cc CX3. Figures [2H1] plot the left- and right-hand sides of (|4.24p - (|4.26p . 
with logx on the horizontal axis. The density Dt was evaluated by numerical 
integration of (I3.8p . using the explicit expressions [HlHI] for 4> and ip. 

Finally, to show the accuracy of the smile asymptotics, we plot the smile 
together with the asymptotic approximations. This is done by simply matching 
Heston prices with Black-Scholes prices by means of a root-finding procedure. 
To evaluate the Heston prices (with initial stock price 5*0 = 1) we use Lee's 
formula [1 21] 



where k is again the log-strike and a is a "damping constant" which we are free 
to choose, noting only that for a > this formula gives us call prices whereas for 
a < — 1 we get the prices of the respective puts. To optimize our results, we will 
use (following Lee) call options for the out-of-the-money strikes, and put options 
for the in-the-money strikes, both with maturity T — I. As a good choice for 
the damping constant a we suggest a = 29.1 for the calls and a — —4.4 for the 
puts. 

The respective Black-Scholes prices are calculated by the Black-Scholes for- 
mula, evaluating the cumulative density function of the normal distribution by 
straightforward numerical integrationo To get good results for deep in-the- 
money/out-of-the- money options, we use as starting point for the root-finding 
procedure the value given by our third order approximation. In the numerical 
example this leads to a stable evaluation of the smile in a quite large interval, 
e.g. log-strikes ranging from —14 to 24. The results, compared with the first- 
and third order asymptotics, are found in Figure [5l There, the log-strike is 
confined to the (more realistic) interval [—2,2]. 



*We thank Roger Lee for helpful comments on this numerical evaluation. 




25 



5000 



10000 15000 20000 25000 



Figure 2: Numerical check for the constant ^3. 
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Figure 3: Numerical check for the constant ^2- 
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Figure 4: Numerical check for the constant Ai. 
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Figure 5: Implied variance cr(fc, 1)^ in terms of log-strikes compared to the first 
order (dashed) and third order (dotted) approximations. 
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